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ABSTRACT 


This report contains useful information for certain aspects of 
the analysis of data obtained by the Ultraviolet Spectrometer and Polarim- 
eter (UVSP) on the Solar Maximum Mission (SMM). The first section describes 
the meaning of the UVSP co-ordinate system and the SMM roll, pitch, and 
yaw and then explains how to overlay UVSP images. The second section 
contains a description of various computer programs that calculate the 
1 ine-of-sight component of the SMM spacecraft velocity from the spacecraft 
to the Sun. The spacecraft velocity can be used to correct or interpret 
the signal observed in UVSP doppl ergrams . This is applied in the third 
section which describes a method of using the spacecraft velocity to 
calibrate UVSP doppl ergrams and magnetograms, i.e., determine the width 
of the observed emission line. The fourth section describes UVSP polariza- 
tion analysis procedures while the fifth section gives expressions for the 
statistical uncertainties in various quantities obtained from UVSP 
measurements. 

Note that the equation numbers in each section begin again at 
Equation (1). 
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CHAPTER 1 

ORIENTING AND OVERLAYING UVSP IMAGES 


1.1 PURPOSE 

The purpose of this chapter is to provide to users of the UVSP 

data explanations of the orientation and coordinate systems applicable to 
UVSP images and of how to overlay different UVSP images and determine their 
position relative to the sun. The UVSP instrument and overall experiment 
have been described by Woodgate et al . (1980) and Henze (1979). 

1.2 DEFINITIONS OF COORDINATES AND PARAMETERS 

The orientation of the SMM and UVSP coordinate systems is shown 
in Figure 1-1. The definition of the UVSP coordinate system is illustrated 
in Figure 1-2. An example of a typical observation showing the signs of 
the various coordinates is given in Figure 1-3 which also shows the defini- 
tion of the angle theta used in the UVSP Data Catalog. 

The pitch and yaw angles and the UVSP X and Y coordinates are given 
in arcseconds. The roll and theta angles are given in degrees. 

SMM Roll Angle - positive when spacecraft is rolled from 

solar north toward west, i.e., clockwise 

SMM Pitch Angle - positive when spacecraft is pointed south 

of sun center when roll angle is zero 

SMM Yaw Angle - positive when spacecraft is pointed east of 

sun center when roll angle is zero 

UVSP X Coordinate - positive toward north when roll angle is zero 

UVSP Y Coordinate - positive toward east when roll angle is zero 

Theta Angle in UVSP Data Catalog - positive from north toward east, 

i.e., counterclockwise. 

The pitch and yaw angles are measured from the center of the solar 
disk and refer to an arbitrary spacecraft "boresight”, the definition of 
which was changed to least twice during the mission. It was, however, 
always located within the lower left quadrant of the UVSP field of view. 

It is possible that up until July 1980 there may have been a center-to- 
limb difference in the location of the spacecraft "boresight". For periods 
of up to one or two days, it is likely that the "boresight" position did 
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Figure 1-1. Orientation of SMM and UVSP Coordinate Systems 
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FIGURE 1-2. Definition and Scale of UVSP Coordinate System 
and Field of View. X and Y are always given in 
integer arcseconds . 
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Figure 1-3. Example of Typical Observation. For this example, 
Pitch is negative and Yaw, Roll, and Theta (e) are 
positive. Of course, X and Y are always either zero 
or positive. 


not change more than a few arcsec which is the probable accuracy of the 
whole SMM pointing scheme. 

The UVSP coordinate system shown in Figure 1-2 is a 256 x 256 arcsec 
array of possible UVSP pointing positions or pixel centers. Because many of 
the entrance slits are larger than 1 arcsec square, the actual possible UVSP 
fiel d-of-view is larger than 256 arcsec square. 

Further quantities associated with a given UVSP data set or experi- 
ment number are the following: 

XCEN or XC = X coordinate associated with the center of a UVSP raster. 

YCEN or YC = Y coordinate associated with the center of a UVSP raster. 

NX = number of X positions in a UVSP raster 
= number of pixels in X direction. 

NY = number of Y positions in a UVSP raster 
= number of pixels in Y direction. 

DX = spacing in arcsec between centers of adjacent pixels in X direction 

DY = spacing in arcsec between centers of adjacent pixels in Y direction 

Note that DX and DY are not necessarily the same as the corresponding 
dimensions of the entrance slit, although they are usually chosen to be the 
same. 

1.3 DATA OBSERVING AND STORAGE ORDER 

In the raster, the value of X or Y corresponding to the outer loop 
would be set at its minimum value for that raster and observations would 
be taken while the value of X or Y corresponding to the inner loop would 
step from its minimum to its maximum value. The outer loop would then step 
to its next higher value and the inner loop would again step from minimum 
to maximum value. The usual nesting of loops had X as the inner loop and 
Y as the outer loop. This meant that the data were taken in vertical strips 
beginning at the lower right. In other words, the very first data point 
which is the same as the first point of the first strip was located at the 
lower right corner of the raster. The last point of the first strip was 
located at the upper right corner. The first point of the last strip was 
located at the lower left corner. The last point of the last strip which 
is the same as the very last data point of the whole raster was located at 
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the upper left corner. In the usual situation where the SMM roll angle 
was zero, north was up and east was to the left. Many display programs 
orient the rasters properly (i.e., the rasters are flipped in the computer). 

1.4 MEANING OF RASTER CENTER COORDINATES 

XCEN and YCEN are stored in Junior's DMA and also in the experiment 
header block of the reformatted data files. The positioning of the raster 
can depend on whether NX, NY, DX, DY are odd or even. The explanation given 
below applies to the X parameters. An identical but separate explanation 
applies to the Y parameters. 

NX Odd: If NX is odd, then the central strip of pixels, i.e., the 

[(NX+l)/2]th strip counting in the X direction, is centered exactly at XCEN. 

NX Even and DX Even: The boundary between the two central strips 

of pixels, i.e., the (NX/2 ) and [(NX/2)+l]th strips counting in the 
X direction, is located exactly at XCEN. The coordinates of the two 
strips are XCEN ± DX/2. 

NX Even and DX Odd: Because the spacing DX is an odd number of 

arcseconds, the real center of the raster - the boundary between the two 
central strips of pixels, i.e., the (NX/2 ) th an d [(NX/2)+l]th strips counting 
in the X direction - is not actually at an integer X coordinate value. The 
boundary between the two central strips is located at XCEN + Thus, the 
X coordinates for the two central strips counting in the X direction are 
XCEN - [ ( DX- 1 )/2] and XCEN + [ (DX+1 )/2] . 

1.5 LOCATION OF EDGE PIXELS IN UVSP COORDINATES 

Because contour plots treat the pixels as points rather than squares 
or rectangles, it is necessary to find the location of the center of the 
edge pixels, i.e., the pointing coordinates for the edge pixels. This 
means that the spread in the coordinates for the edge pixels at opposite 
sides of the raster is not the same as the field of view because the field 
of view includes the size of the entrance slit. 
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X Coordinate for Edge Pixels 



Upper 

Lower 

NX Odd 

XCEN + * ox 

XCEN - * DX 

NX Even, 
DX Even 

XCEN + (NX-1) * y- 

XCEN - (NX-1) * y 

NX Even, 
DX Odd 

XCEN - ( D - X ~-)- +(y- * Dx) 
or 

XCEN + 7f [(NX-1 )*DX+1] 

XCEN + * Dxj 

or 

XCEN - ( DX ~ 1 ) -(^r-i)* ox 
or 

XCEN - ^ [(NX-1) * DX - 1] 


Again, a similar set of expressions involving NY and DY is applicable for 
the Y coordinate of the left and right edge pixels. 

1.6 OVERLAYING RASTERS IN UVSP COORDINATES 

When feature tracking was being used and the rasters were not too 
far separated in time, one can simply determine the X and Y coordinates 
of the upper left corner pixel in each raster, subtract the underlay 
coordinates from the overlay coordinates to find the position of the 
upper left pixel of the overlay relative to the upper left corner of the 
underlay, and then place the overlay on the underlay using the position 
in arcsec of the upper left corner of the overlay relative to the upper 
left corner of the underlay. Of course, the scales of both rasters should 
be the same. To determine the scale, one should recall that the number of 
arcsec between the pixel centers on the top edge and on the bottom edge 
is (NX-1)*DX; a similar expression applies in the Y direction. On the 
other hand, some contour plots are surrounded by a box that corresponds 
to dimensions of NX*DX and NY*DY; in such plots the contour lines do not 
fill up the entire box, i.e., there are blank strips one pixel wide along 
the adjacent edges. 


6 










1.7 


OVERLAYING RASTERS IN PIXEL NUMBERS 


Because contour plots are sometimes labelled in pixel numbers, 
it is then convenient to use the underlay pixel number scale as the 
coordinate system. One goes through the process described in Section 1.6 
above and then divides the offsets by the underlay values of DX and DY. 

If the pixel number scale has its origin at the upper left and if the 
contour plots are pushed into the upper left corner (i.e., if the blank 
strips are at the bottom and right edges), then the upper left corner 
of the overlay can be positioned directly on the underlay using the 
underlay pixel number scale. (Of course, the upper left corner of the 
overlay could be outside the box on the underlay.) 

1 .8 CORRECTING FOR SOLAR ROTATION: SPACECRAFT PITCH AND YAW 

If feature tracking was not being used, it is necessary to correct 
for solar rotation. One should calculate the number of arcsec by which a 
feature would have shifted to the west (usually to the right) in the later 
raster. If the pitch and yaw were unchanged, then one would simply shift 
the later raster to the east (usually left) by the appropriate number of 
arcsec caused by rotation. If the pitch and yaw were changed, then one 
should also shift one of the rasters in the proper direction. 
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CHAPTER 2 

SPACECRAFT VELOCITY PROGRAMS - COMPONENT 
OF ORBITAL VELOCITY ALONG LINE OF SIGHT TO SUN 

2.1 INTRODUCTION 


This chapter describes Fortran and IDL programs and subprograms 
and associated data files which are utilized to calculate the component of 
the orbital spacecraft velocity along the line of sight to the sun at any 
time during the useful part of the SMM mission in 1980. 

Definitions for some of the symbols are as follows: 

VLOS = Line-of-sight component of velocity relative to sun 

in km s _1 . Negative when spacecraft is moving toward 
the sun, i.e., before orbital noon; positive when moving 
away from the sun, i.e., after orbital noon. Floating 
point. 

IDOY = Day of year for which velocity is desired. Integer. 

FRD = Time expressed as fraction of day for which velocity is 
desired. Floating point. 

IOPEN = Code in some versions of the IDL procedures. It indicates 
whether the call is the first one (IOPEN = 0) or a second 
or following call (IOPEN f 0). During the first call to 
VELSCI, IOPEN is set equal to a non-zero number so that 
if the calling program or procedure uses a variable 
instead of a constant for that parameter, then the same 
variable can be used in following calls without changing 
its value; it must have been set equal to zero before 
the first call. Integer. 

Someone who wants to use the Fortran subroutines and IDL procedures 
without learning any background can read Sections 2.2 and 2.3 only. The 
explanation of the method is given in Section 2.4. Section 2.5 describes 
the programs used to create the data files containing the tables of orbital 
information required; Section 2.6 describes the data files themselves. 
Section 2.7 describes the different versions of the Fortran subroutines 
and the IDL procedures available for the calculation of the line-of-sight 
velocity relative to the sun. 
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2.2 


FASTER AVAILABLE VERSIONS OF THE VELOCITY CALCULATION SUBPROGRAMS 
FOR INDIVIDUAL TIMES - SHORT LISTING 


The subprograms listed here open and read in the data file on the 
first call only. The second and following calls should be much faster. 

For explanation, see the later sections. 

2.2.1 Fortran 

VELSC2.FTN - contains SUBROUTINE VELSC(IDOY ,FRD,VL0S). 

This version of VELSC uses data file BETAT1.DAT. 

2.2.2 IDL 

It is not known which of the versions in the two files given here 
is faster. All parameters must be scalars. 

VELSC2.PR0 - contains PRO VELSCI ,IDOY,FRD,VLOS,IOPEN. 

This version of VELSCI uses data file BETAT1.DAT. 

VELSC4. PRO - contains PRO VELSCI ,IDOY,FRD,VLOS,IOPEN. 

This version of VELSCI uses data file BETAT3.DAT. 

2.3 IDL VELOCITY CALCULATION PROCEDURES TO HANDLE ONE-DIMENSIONAL ARRAYS 
FOR INPUT PARAMETERS - SHORT LISTING 

It is not known which of the versions given here is faster. The 
time FRD (in fraction-of-day units) and the velocity VLOS can be one-dimen- 
sional arrays (same size); the times in FRD sholild be all in the same orbit. 
VLOS is always a scalar. 

VELSCI. PRO - contains PRO VELSC, IDOY, FRD, VLOS. This version 
of VELSC uses data file BETAT1.DAT. 

VELSC3.PR0 - contains PRO VELSC, IDOY, FRD, VLOS. This version 
of VELSC uses data file BETAT3.DAT. 


2.4 METHOD 

It is assumed that the spacecraft orbit is circular and that the 
orbit of the earth is circular. We thus ignore any motion of the earth 
toward or away from the sun. At the present time, the circular spacecraft 
velocity is assumed to be constant over the entire mission. 
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The equation for the component of the spacecraft velocity along 
the line of sight to the sun, V||, is 

V || = V 0 sin £ 

where t Q = time of spacecraft noon for this orbit, 
t = time for which V|| is being calculated, 

T = orbital period of spacecraft, 

V 0 = spacecraft orbital velocity, 

3 = angle between spacecraft orbital plane and direction toward 
sun. 

The time of spacecraft noon is interpolated from a table with entries 
usually spaced every 7 to 10 orbits. The value of 3 is also interpolated 
from the same table. The table is based on the orbital predictions contained 
in the blue books in the EOF at GSFC. The relevant information is on the 
right hand side of the pages of the Equator Crossing List which usually 
appears in both Volumes I and II of the blue books. Orbital noon is the 
same as the time of maximum solar elevation. 


j cosb 


2.5 PROGRAMS TO GENERATE DATA TABLES 
VSCDAT1 . FTN, .TSK 

Fortran program to create and/or add to file BETAT1.DAT; 
inserts new data from keyboard into proper location in 
table in file. Old version of file BETAT1.DAT should be 
purged after new data have been inserted. The FTN file 
contains the statement INCLUDE DOY.FTN. 

DOY.FTN 

Fortran SUBROUTINE D0Y(ID0Y, IYR, IMON, I DOM) . Calculates 
day-of-year for input year, month, day-of-month. Used by 
VSCDAT1 because blue book with times of orbital noon and 
8 gives month and day-of-month instead of day-of-year. 

VSCDAT3 .PRO. 

IDL program which creates file BETAT3.DAT based on input 
file BETAT1.DAT to allow IDL ASSOC function to be used for 
reading in data table. This program should be run every time 
file BETAT1.DAT is changed and the old version of BETAT3.DAT 
should be deleted before running VSCDAT3.PR0 or purged 
afterward. 
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2.6 FILES CONTAINING DATA TABLES 


BETATI.DAT 

ASCII data table containing 6 angle and time of orbital noon 
for whichever orbits were entered. For 1980 only. Order of 
entries is according to time. Each record contains p 
(in degrees), day-of-year, hr, min, sec with format (FI 0.3, 
417). Can be printed out with PIP. At the present time, 
the spacing between entries is usually 7 to 10 orbits; 
however, additional orbits can always be inserted using 
VSCDAT1 . An editor can also be used to insert additional 
entries or delete mistakes; however, care must be taken to 
insert the new entry into the proper position according 
to the time. 

BETAT3.DAT 

Data table created by IDL program VSCDAT3.PR0 and containing 
primarily p angle and time of orbital noon for whichever 
orbits were in BETATI.DAT. Each record contains three 
unformatted (binary) floating point numbers. Each number 
in the first record (IDL record 0) is equal to the number 
of entries of p and time in the table. Let the number of 
entries be NSIZ. Then the total number of records in 
BETAT3.DAT = NSIZ + 1 which is one greater than the number 
of records in BETATI.DAT. For the second and following 
records in BETAT3.DAT (IDL records 1 through NSIZ), the 
first number is P in degrees, the second number is the 
day-of-year (remember, in floating point), and the third 
number is the time in fraction-of-a-day . The reason for 
storing the DOY and fraction-of-day separately is to preserve 
significant figures. This file can be read with the IDL 
ASSOC function. 


2.7 SUBPROGRAMS WHICH CALCULATE THE LINE-OF-SIGHT COMPONENT OF 
SPACECRAFT VELOCITY RELATIVE TO THE SUN 

The Fortran subroutines and IDL procedures all calculate the line- 
of-sight velocity component VLOS (in km s -1 ) for a given day-of-year IDOY 
and fractional time-of-day FRD in 1980. 

2.7.1 Fortran Subroutines 

Both files contain a SUBROUTINE VELSC( IDOY, FRD, VLOS) . 

Both subroutines use data file BETATI.DAT. 


11 


VELSCT.FTN 

For each call to VELSC, opens file BETAT1.DAT and reads 
individaul records until the required records are found. 

Slow. 

VELSC2.FTN 

On first call only to VELSC, opens and reads in entire file 
BETAT1.DAT which is then saved in memory. Search through 
table on following calls (and also on first call) is made 
by searching through table already stored in memory. Faster 
for second and following calls. 

2.7.2 IDL Procedures 

The files contain either PRO VELSC, IDOY ,FRD,VL0S or PRO VELSCI, 
ID0Y,FRD,VL0S,I0PEN which do same thing as Fortran Subroutine described 
above. Procedure VELSCI includes the parameter IOPEN to indicate whether 
the call is the first one (IOPEN = 0) or a second or following call 
(IOPEN f 0); on the first call, the data file is opened and read in its 
entirety, and the resulting table is saved in COMMON BLOCK so that it 
does not need to be read in again on following calls. The parameter IOPEN 
is set equal to a non-zero number in the procedure VELSCI during the first 
call so that if the calling program or procedure uses a variable for that 
parameter instead of a constant, then the same variable can be used in the 
following calls without changing its value; of course it must be set equal 
to zero before the first call. Procedure VELSC opens data file and reads 
individual records for each call and should therefore be slower than 
VELSCI for the second and following calls for individual times. However 
VELSC can handle one-dimensional arrays (same size) for FRD and VLOS; 
the times in FRD should be in the same orbit. The velocity is calculated 
using IDL array arithmetic which should be faster for arrays. IDOY is 
always a scalar. 

VELSCI. PRO - contains VELSC (IDL procedure). Uses BETAT1.DAT. 

Slow for individual times but can handle one-dimensional arrays 

for FRD and VLSO. 

VELSC2.PR0 - contains VELSCI (IDL procedure) . Uses BETAT1.DAT. 

Faster for second and following calls. FRD and VLOS are scalar. 

VELSC3.PR0 - contains VELSC (IDL procedure). Uses BETAT3.DAT. 

Slow for individual times but can handle one-dimensional arrays 

for FRD and VLOS. 

VELSC4.PR0 - contains VELSCI (IDL procedure). Uses BETAT3.DAT. 

Faster for second and following calls. FRD and VLOS are scalar. 
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2.7.3 SOL Procedure 


One of the IDL procedures was converted to the language SOL. 

The procedure is in the file VELSC1 SOL .PRO and is effectively identical 
to the procedure in the IDL file VELSC1.PR0. All of the input and 
output parameters, etc., are identical. 
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CHAPTER 3 


ANALYSIS OF UVSP DOPPLERGRAMS AND MAGNETOGRAMS AND THEIR 
CALIBRATION USING THE ORBITAL VELOCITY OF THE SMM SPACECRAFT 

3.1 INTRODUCTION 

The provision in the Ultraviolet Spectrometer and Polarimeter (UVSP) 
experiment (Woodgate et al . 1980) on the Solar Maximum Mission (SMM) space- 
craft of pairs of exit slits allowed different parts of one emission line 
(red and blue wings) to be observed simultaneously. Comparison of the 
intensities in the two parts of the line made it possible to observe 
relatively rapidly the 1 ine-of-sight components of the velocity (the 
dopplergram observing mode) and, when the polarimeter was also being used 
to measure the circular polarization caused by the Zeeman effect, the 
magnetic field (the magnetogram mode). The magnetogram mode also produced 
the same velocity information as the dopplergram mode. A procedure for 
analyzing dopplergrams and magnetograms obtained with wide exit slits has 
been described by Henze et al . (1982). 

To convert the observed dopplergram and magnetogram signals to 
line-of-sight components of the velocity and magnetic field, one must know 
the width of the emission line for each pixel. (We assume a gaussian line 
profile where the width is specified by the doppler width. We also assume 
that there is no continuum emission.) A line-width calibration feature was 
included in the control system for the instrument. This feature, which was 
not always activated during the period of SMM operations in 1980, allowed 
periodic shifts by a known amount of the spectrometer grating and, therefore, 
of the line relative to the exit slits. The known shift thus provided an 
artificial velocity to calibrate the measured signal. A method of analysis 
using this feature was also described by Henze et al . (1982). When the 
feature was not used, it is still possible under certain conditions to 
determine the line width using the line-of-sight component of the orbital 
velocity of the spacecraft as the calibration reference. The primary 
purpose of this paper is to describe this method which uses the technique 
of least squares to fit the observed doppler shift signal for a pixel to the 
known line-of-sight component of the spacecraft velocity. The fitting pro- 
cedure yields the doppler width and the offset for the pixel of the mean 
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position of the line relative to the position of the exit slits. The 
method is feasible because the maximum doppler shift due to the spacecraft 
motion is a significant fraction of the typical solar line widths. 

There are several necessary conditions for application of the 
method. The observational data should consist of many repetitions of the 
dopplergram or magnetogram raster with the total duration covering a large 
fraction of the spacecraft orbit so that the 1 ine-of-sight component of the 
spacecraft velocity varied significantly during the observing sequence. 

There should have been no large secular or long-duration changes in the 
velocity of the solar emitting region during the observing sequence. 

Finally, the line width must have been constant during the observing 
sequence. 

As background for the description of the calibration method and to 
show where the doppler width is required, we also summarize the theoretical 
expressions for conversion of the observed signals to velocities and 
magnetic fields. Although, as mentioned above, the expressions for wide 
exit slits have been published previously, those for narrow exit slits 
have not. 

3.2 THEORETICAL INTERPRETATION OF DOPPLERGRAM AND MAGNETOGRAM SIGNALS 

There are two types of exit slits used for dopplergrams and magneto- 
grams in the UVSP: wide and narrow slits. Each member of a wide pair is 

0.3 A wide in second order; the two members are exactly adjacent so that 
the entire line is effectively covered by the pair of slits. Each member 
of a narrow pair is 0.03 A wide in second order; the centers of the two 
members are 0.294 A apart at a wavelength of 1548 A (Athay et al . 1982). 

An illustration of the poition of the wide and narrow exit slits relative 
to line profiles with two different widths is given in Figure 1 of the 
paper by Gurman and Athay (1983). In the remainder of this paper, the 
wide slits will be assumed to be semi-infinitely wide (so that together, 
the pair always cover the entire line) and the narrow slits assumed to be 
infinitesimally wide. 
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In the UVSP dopplergram or magnetogram observation, the observed 
doppler shift signal is defined as 


R v = 


*r ~ *b 
Ir + lb ’ 


0 ) 


where I is the intensity (integrated across the exit slit) and r and b 
denote the exit slits for the red and blue wings of the line. In a UVSP 
magnetogram observation, the observed magnetic circular polarization 
signal is defined as 





( 2 ) 


where V is the Stokes parameter for circular polarization. 

The 1 ine-of-sight velocity v is related to the line shift A A v by 


v = (c/A 0 )AA v . 


(3) 


In addition to this doppler shift, any offset due to the non-ideal 
positioning of the line relative to the pair of exit slits will add a 
constant offset term to the line shift and thus to the velocity. For 
wide slits (Tandberg-Hanssen et al . 1981), 

AA V = AAq inverf(Ry) > (4) 

where AAq is the doppler width and inverf is the inverse error function. 
If the shift is small (AA v /AAd « 1), 

AAy ~ AAq(u 2 /2) Ry . (5) 

For narrow si i ts , 


AA V = (AAp/AA sep ) arctanh(R v ) , (6) 

where AA se p is the wavelength separation between the two exit slits and 
arctanh is the inverse hyperbolic tangent function. If the shift is small, 

AAy ~ (AAg/AAggp) Ry . (7) 
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The magnetic field B is related to the Zeeman splitting AAg by 


n . 4ir m c 2 
B - ax B • 


For wide slits (Henze et al . 1982), 

AAg cosy = AA D (iri/2) R m [1 - erf 2 (AX v /AXQ)] exp(AA v /AA D )‘ 


( 8 ) 


(9) 


where y is the angle between the field and the line of sight. If the line 
shift is small. 


AAg cosy - AXq(tt^/ 2) R m . (10) 

For narrow slits, 

AAg cosy = (AXp/AX se p) R m , (11) 

which is independent of the line shift. 

3.3 CALIBRATION USING SPACECRAFT ORBITAL MOTION 

The basic method consists of fitting the appropriate function 
(inverf or arctanh) of the observed values of the doppler shift signal R v 
for a pixel to a linear function of v sc , the time-varying component of the 
spacecraft orbital velocity along the line of sight from the spacecraft 
to the sun. The spacecraft velocity component is calculated from the 
orbital ephemeris for the time of each observation as explained in 
Chapter 2 of this report. Maximum values are approximately ±7.5 km s -1 . 
For wide slits, 

inverf (R v ) = a 0 + ai v sc (12) 

For narrow slits, 

arctanh (R v ) = a 0 + ai v sc . (13) 
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The least squares fitting procedure yields the parameters a 0 and 
a t ; their uncertainties are proportional to the standard deviation of the 
difference between the observed function of R v and the fitted curve 
(Bevington 1969). Then, for wide slits, 

AXq = W(cai) , (14) 


and, for narrow slits, 

AAq = [ ( A 0 AA se p)/ca i ) ]^ . (15) 

For both types of slits, 


v Qff = a 0 /ai , (16) 

where the offset velocity v 0 ff includes both any steady 1 ine-of-sight 
velocity on the sun and the component due to the possible non-ideal position 
of the exit slits relative to the line. The uncertainties for the doppler 
width and the offset velocity can be obtained by propagating the uncertain- 
ties for a 0 and aj through Equations 14 to 16. 

The method described above has been incorporated in an IDL procedure 
for convenient calculation on a computer. 

3.4 APPLICATION TO WIDE-SLIT MAGNETOGRAMS 

As a demonstration of the method, it has been applied to two 
magnetogram observations analyzed previously by Henze et al . (1982). The 
data were obtained in the C IV line at 1548 A in the umbra of a sunspot on 
23 and 25 October 1980; the exact locations are shown in Figure 1 of the 
paper by Henze et al . 

Figure 3-1 shows the fit for one pixel in one of the magnetograms 
(UVSP Experiment 14398) of the observed doppler shift signal (converted 
to velocity) to the 1 ine-of-sight component of the spacecraft velocity. 

The complete sequence displayed here consisted of 256 repetitions of 
5.056 s each. The observed curve shows fluctuations due to statistical 
noise (high frequency, small amplitude) and umbral oscillations (-150 s 
period, -3 km s -1 amplitude). Although both sources of fluctuations 
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contribute to the formal uncertainties in the resulting parameters, the 
good fit of the general trend of the observed curve to the spacecraft 
curve and the large number of measurements causes those uncertainties 
to be smal 1 . 

The numerical results for all five pixels are given in Table I. 

The doppler widths and the offset velocities along with their uncertainties 
(1 sigma) were obtained as described above. The offset velocities are all 
large and of the same sign, indicating that the exit slits were not well 
centered on the line. The measured values of the doppler width were then 
used to infer improved values of the magnetic field from the observed 
circular polarization signals; previously, the widths had been assumed 
to be 100 mA for all of the pixels. Because the values of the inferred 
fields are not greatly changed, the conclusions of Henze et al . (1982) 
about the existence of fields of the order of 1000 gauss in the transition 
region above a sunspot are not affected. 

Table I also compares two different approaches to determining the 
magnetic field. In one approach, the full expression in Equation 9 was 
used for calculating the Zeeman splitting, and the magnetic field was 
determined separately for each repetition. The mean value is given in 
Table I where the uncertainty is the standard deviation of the mean. 

Using the other approach, Henze et al . (1982) summed the observed count 
rates over all 256 repetitions at each pixel and waveplate position before 
performing the magnetic field analyses; they also used the approximation 
in Equation 10. The uncertainty was derived from the statistical 
uncertainty in the count rates. Their results are also given in Table I 
after correction for the doppler width obtained here. The two methods 
yield similar values for the magnetic fields which agree well within the 
uncertainties. However, it should be better in principle to use the first 
approach, determining the magnetic field separately for each repetition 
and then taking the mean. The uncertainties also agree closely, thus 
indicating that the main source of random error in the magnetic field 
determination in this data set is the statistical counting noise and not 
any real intensity or magnetic field fluctuations on the sun. Of course, 
in other observations where the count rates are higher and the statistical 
noise in the results lower, real fluctuations may contribute significantly 
to the uncertainties. 
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TABLE I 


Offset Velocity, Doppler Width, and Longitudinal 
Magnetic Field, UVSP Experiment 14398 


Pixel 

v off 
(km s -1 ) 

(mA) 


B cos y* 
(gauss) 

B COS Y** 
(gauss) 

1 

16.4 

± 0.7 

94 

+ 

4 

-1320 

+ 

630 

-1090 ± 620 

2 

14.4 

± 0.6 

88 

+ 

4 

-320 

± 

500 

-370 ± 480 

3 

10.1 

± 0.5 

80 

± 

3 

-1010 

± 

430 

-830 ± 360 

4 

10.6 

± 0.5 

84 

± 

3 

-960 

± 

330 

-940 ± 290 

5 

11.8 

± 0.6 

89 

± 

4 

-1060 

± 

320 

-1010 ± 300 


* Mean of results using Equation 9 for each repetition. 

** From Henze et al . (1982) but corrected for doppler width 
given here instead of the previously assumed value of 100 m&. 
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As a check on the method, it was also applied to a calibrated 
magnetogram (UVSP Experiment 14302). The sequence contained 135 repetitions 
of 15.5 s each. Table II compares the results for those pixels that showed 
no systematic changes during the observations except for periodic oscilla- 
tions. The widths obtained by Henze et al . (1982) were from the analysis 
of the first occurrence of the special calibration shift in the sequence 
(the 16th repetition). The results agree well. The small uncertainties 
in the results from this paper are due to the inclusion of many more 
repetitions in the analysis. 


TABLE II 

Comparison of Methods of Calibration of the Line Width, 
UVSP Experiment 14302 


ixel 

AAp (mA) 


Pixel 

AAg 

(mA) 


Pixel 

AA d 

M) 



★ 

** 


* 




•k 


★★ 

1 

144 ± 3 

135 

8 

134 ± 

2 

139 

15 

136 ± 

2 

122 

2 

142 ± 2 

137 

9 

152 ± 

2 

150 

16 

137 ± 

1 

126 

3 

109 ± 2 

111 

10 

109 ± 

1 

100 

17 

112 ± 

1 

104 

4 

116 ± 3 

124 

11 

101 ± 

2 

100 

18 

121 ± 

1 

117 

* From 

method des 

cribed 

in this 

paper. 








** From Henze et al . (1982), based on special calibration feature. 
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3.5 COMPUTER PROCEDURES 

Some of the calculations necessary to use the above method have 
been incorporated in procedures written in the languages IDL and SOL. 

The procedures are in the files: 

IDL - LSVCALIDL.PRO 
SOL - LSVCALSOL.PRO. 

Both files contain the procedure: 

PRO LSVCAL , VSC , JSL IT , FRV , CV , DL D , VOFF , S I 6DLD ,SIGVOF,SIG. 

Also needed are the procedures: 

PRO INVERF, X, Y and 
PRO ARCTNH, X, Y. 

They are included in the SOL file named above and are automatically 
compiled when the IDL file is compiled. The IDL versions are in. the files 
INVERF. PRO and ARCTNH. PRO. Comments included in the files explain all of 
the parameters. 

3.6 CONCLUSIONS 

A method of determining the doppler width for each pixel and 
thereby calibrating UVSP dopplergrams and magnetograms using the line-of- 
sight component of the spacecraft orbital velocity has been described and 
illustrated. The method is useful when the special UVSP calibration feature 
was not activated, as was often the case during the period of SMM operations 
in 1980. The method requires the assumptions of constant doppler width and 
constant mean doppler shift (i.e., 1 ine-of-sight component of the flow 
velocity) during the observing sequence. The sequence must be long enough 
to determine the parameters to the desired degree of accuracy. If there 
are any oscillations present, then the duration should include several 
complete periods. Because of these restrictions, a better approach in 
the event that future SMM/UVSP observations become possible would be to 
invoke the normal calibration shift on every second raster, thus assuring 
an almost continuous calibration for the line widths. This procedure was 
used occasionally during 1980 and was called the Double Dopplergram 
Determination by Simon et al . (1982). 
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CHAPTER 4 

UVSP POLARIZATION ANALYSIS PROCEDURES 


This cahpter describes the analysis of UVSP polarization obser- 
vations. Such observations are required, for example, to measure the 
longitudinal magnetic field as described in another chapter. 

4.1 POLARIZATION SIGNAL 

For a polarized light source described by the Stokes parameters, 

I, Q> U, V (also sometimes denoted by S 0 , Si, S 2 , S 3 ), the rotating 
waveplate in the UVSP modulates the detected signal $ as a function of 
the angle to of the waveplate fast axis: 

$(to) = K[I + p(a cos28 + b sin20 sin4io + b cos20 cos4to)Q + 

+ p ( -a sin20 - b cos20 sin4to + b sin20 cos4to)U + 

+ p sin6(cos20 sin2to - sin20 cos2to)V ], (1) 

where K is the instrumental efficiency for the detector (including detector 
sensitivity, mirror and grating reflectivities, telescope collecting area, 
entrance slit solid angle, integration time, and waveplate transmissivity), 
p the degree of polarization of the analyzer (grating in this case), 0 the 
angle of analyzer axis, 6 the retardance of the waveplate, a = (1 +cosS)/2, 
and b = (1 -cos6)/2. The angles w and 8 are both positive clockwise (the 
direction of the waveplate rotation) when looking toward the light source 
and are measured relative to the same arbitrary reference direction, 
chosen to be the positive UVSP X axis in the instrument so that 0=0. 

In the usual spacecraft roll orientation, the positive X axis pointed to 
solar north. At the wavelength of 1548 A, p = -0.65 and 6 = 270°. The 
waveplate position at which the first data point is measured in any UVSP 
polarization observation is at w = 315°. Equation (1) is a Fourier series 
and the observed values of $(w) can be Fourier analyzed to determine the 
coefficients and thus the Stokes parameters. 
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4.2 


MAGNETOGRAMS 


In magnetogram observations, only four equally spaced measurements 
were made so that the 4w components (and thus the linear Stokes parameters, 
Q and U , which are expected to be small with the wide exit slits, cannot 
be determined. Also, the constant (Oto) paQ term is ignored because it is 
expected that Q « I . Equation (1) then simplifies to 

$(oj) = K[I + p sirnS V sin2oi] . (2) 

The Fourier analysis yields 

KI = [$(315°) + $(45°) + $(135°) + $(225°)]/4 (3) 

and 

KV = -[$(315°) - $(45°) + $(135°) - $(225°)]/(4p sinS). (4) 
The fractional circular polarization is then V/I. 

4.3 POLARGRAMS 

Although UVSP polargrams can be made with any number of waveplate 
positions, the most common experiment definition would have sixteen equally 
spaced measurements (i.e., 22.5° apart). For 9 =0, Equation (1) becomes 

$(oo) = K[I + paQ + psinfi V sin2w 

+ pbQ cos4u> - pbllsin4w] . (5) 

To make the following equations easier to understand, let us denote 
the measurements at different waveplate positions as follows: 


$1 

= $(315°), 

$2 

= $(337.5°), 

$3 

= $( 0° ), 

4>4 

= $(22.5°), 

$5 

= $(45°), 

$6 

= $(67.5°), 

$7 

= $(90°), 

$ 8 

= $(112.5°), 

4*9 

= $(135°), 

$10 

= $(157.5°), 

*11 

= $(180°), 

$12 

= $(202.5°), 

$13 

= $(225°), 

$14 

= $(247.5°), 

$15 

= $(270°), 

$16 

= $(292.5°). 
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The Fourier analysis yields 



16 

E $i 
i=l 


+ a[<X> 1 - $3 + $ 5 - $7 + $9 - $ u 

+ $13 ~ $15]/(8b) , 

which, if the polarization is small, simplifies to 

1 16 

KI ~ tV E $i . 

16 i = l 


The Fourier analysis also yields 

KQ = [$! - $3 + $ 5 - $7 + $9 - $ U + $ 13 - $is]/(8pb), 
KU = f$2 “ $4 "*■ $6 ~ $8 $10 “ $12 $14 “ $16 ]/(8pb) , 


and 


KV 


- j [$i - $ 5 + $9 - $13 3 

+ [ $2 ~ $4 “ $6 + $8 

+ $io - $12 - $14 + $16 ]/ 2 * |/ 8 p sin 6 ) . 


The fractional linear polarization is then 

P = (Q 2 + U 2 )*/I 

and the angle V of the linearly polarized component relative to the 
X-axis is given by 

tan IV = U/Q . 


( 6 ) 

( 7 ) 

( 8 ) 
( 9 ) 

( 10 ) 

( 11 ) 

UVSP 

( 12 ) 
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CHAPTER 5 

STATISTICAL UNCERTAINTIES 


This chapter gives expressions for the statistical uncertainties 
in UVSP measurements of various quantities. The quantities themselves 
were mostly defined earlier in Chapter 3. The only uncertainty considered 
is that due to the fluctuations inherent in measuring a number of events 
such as photons counted by a detector. 


5.1 BACKGROUND 

The basis for the determination of an uncertainty is that the 
uncertainty e in the measurement of a number of counts $ is given by 

e($) = . (1 ) 

The uncertainty used here is, therefore, the standard deviation (also often 
referred to as one sigma). 

The uncertainty in the sum or difference of two quantities whose 
uncertainties or errors are independent is given by 

e(A + B) = [e 2 (A) + e 2 (B)]± (2) 

and 

e(A - B) = [e z (A) + e 2 (B)]* . (3) 

The uncertainty in the sum or difference of two quantities whose 
uncertainties or errors are not independent is given by 

e(A + B) = e(A) + e(B) (4) 

and 

c(A - B) = e(A) + e(B) . (5) 


Similar expressions are valid for the relative uncertainties of 
products and quotients. For independent uncertainties, 


and 


e(AB) = | AB 


e 2 (A) , c 2 (B) 
A 2 B 2 


i 


e(A/B) = 


| A | e 2 (A) 

i B l A 2 



( 6 ) 

(7) 
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( 8 ) 


For uncertainties or errors that are not independent, 



£(AB) = | AB | 

e(A) 



. I A I 

l B l . 

and 





e(A/B> = |f| 

e(A) 

£ ( B ) 1 


. |A| 

|B| . 

A special case is 





£ ( A 2 ) = 2 | A | £ 

(A) ; 

the inverse is 





:(A±) = . 

2A* 


If a quantity is multiplied by a constant a, the uncertainty is 

e(aA) = a e (A) . 


(9) 


( 10 ) 

(ID 


( 12 ) 


5.2 POLARIZATION 

The uncertainty in the intensity Stokes parameter depends on how 
many waveplate positions were used in the measurement. (Polarization effects 
are assumed to be small.) If N is the number of waveplate positions, 

e(KI) = 4 (SU* » 03) 

N* 

where 

N = 1 for dopplergrams, 

N = 4 for most magnetograms, 

N = 16 for most polargrams. 

If N = 4, 

e(KV) = £ ( KI ) / ( p sin6) . (14) 


If N = 16, 


e(KV) = 2* e(KI)/(p sin6) , 

(15) 

e(KQ) = 2i £ ( KI ) / ( pb ) , 

(16) 

e(KU) = e(KQ) . 

(17) 


and 
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Also, 


e(K[Q 2 + .U 2 ]*) = e(KQ) , 


(18) 


and e(H') = i arc sin 

If the linear polarization is small. 


e(KQ) 

K[Q 2 +U 2 ]i • 


e(P) = e([Q 2 + U 2 ]i)/I = e(KQ)/(KI) . 


(19) 


( 20 ) 


5.3 LINE INTENSITIES 

The combined line intensity I|_ is the sum of the two intensities 
I r and I b : 

l L = l r + l b , (21) 

where r and b denote the detectors in the red and blue wings, respectively. 
(For wide exit slits, this is the true total line intensity.) Because we 
are actually measuring counts (e.g., K r I r instead of I r ), we define the 
weighted sum of the two count rates: 

[K r K b ]± II = R _i K r I r + Ri Kblb » (22) 

where 

R = K r /K b . (23) 

. The uncertainty is then 

c([K r K b ] i I L ) = [R- 1 e 2 (K r I r ) + R e 2 (K b I b )]*-. (24) 


5.4 MAGNETOMETRY 

If the polarization is small, 

c(V/I) = c(V)/I = e(KV)/(KI) . (25) 

The uncertainty in the observed magnetic circular polarization signal is 


e(Rrn) " 2 


c 2 (KrVr) + e 2 (K b V b ) 


(K r I r ) 2 


(K b I b ) 2 J 


(26) 


For N = 4, 


£(Rm) 


1 + 1 


K r I r I b 


/ (4 p sin 6) 


(27) 
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For N = 16, 


e(R m ) = 2i 


Kr^r Kb I b 


/ (8 p sinS) . 


(28) 


Allowing for any uncertainty in the doppler shift (such as that resulting 
from the least-squares fitting technique using the spacecraft velocity 
described in Chapter 3) and assuming that the doppler line shift is small, 
we obtain, for wide exit slits. 


c(AAg cosy) = |AAg cosy 


e 2 ( AX p) 

. a *d 2 


+ 



For narrow exit slits. 


e(AAg cosy) = |AAg cosy 


4e 2 (AAp) + e 2 (R m ) ' 


aa d 2 


R m 


For either type of exit slit. 


e(B cosY) 
B cosy | 


e(AAg cosY) 
| AAg cosy | 


(29) 


(30) 


(31) 


5.5 VELOCITY MEASUREMENTS 

If the uncertainties in the numerator and denominator of the observed 
doppler shift signal are not assumed to be small and if the doppler shift 
itself not assumed to be small, the uncertainty in the shift signal is 


e(R v ) = (l + |R v l) 


e([K r K b ]* I L ) 
[K r K b ]^ I b 


(32) 


Including the uncertainty in the doppler width yields, for wide exit slits, 

li 


e(AA v ) = j AA V 


e 2 ( AA p) + £ 2 (R V ) 

AA 0 2 R v 2 


(33) 
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For narrow exit slits, 


£ ( AAy ) - AAy | 


' 4 £ 2 (AAp) + 
AA D 2 



For either type of exit slit, 

e (v) = £(AAy) 


(34) 


(35) 
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